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ABSTRACT 

•  An  interdisciplinary  research  group  in  modeling,  analysis  and  simulation  was  - 
established  at  Florida  A&M  University  (FAMU).This  group  set  as  its  scientific 
goal  the  development  of  robust  computational  algorithms  that  £ire  weU  suited  for 
nonconvex  variational  principles,  and  other  nonlinear  problems  that  arise  in  the 
study  of  highly  nonlinear  materials.  The  equally  important  second  goal  of  the 

•  group  was  to  enhance  the  participation  of  under-represented  minorities  in  applied 
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•  relationship  between  FAMU,  Carnegie  Mellon  University,  University  of  Minnesota, 
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research  goeds  of  this  project,  and  delineates  the  accomplishments  made  thus  far. 
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1.  Introduction 


Florida  A&M  University  has  established  an  interdisciplinary  group  in  modeling, 
analysis,  and  simiUation.  The  main  scientific  goal  of  this  group  is  the  develop¬ 
ment  of  robust  computational  algorithms  for  nonconvex  variational  principles, 
and  other  nonlinear  problems  that  arise  in  the  study  of  highly  nonlinear  materi¬ 
als.  This  group  was  composed  of  James  C.  Turner  Jr.  and  Yanzhao  Cao  firom  the 
Department  of  Mathematics  at  Florida  A&M  University,  and  consultants  Richard 
James  (University  of  Minnesota),  David  Kinderlehrer  (Carnegie  Mellon  Univer¬ 
sity),  Mitchell  Luskin  (University  of  Minnesota),  and  Gregory  Forest  (University 
of  North  Carolina). 

An  additional  goal  for  Florida  A&M  University  (FAMU)  was  the  establishment 
of  a  Center  for  Computational  Science  (CCS).  The  aim  of  this  Center  was  to  pro¬ 
vide  world-class  education  and  research  in  areas  that  are  essential  to  the  nation’s 
scientific  progress  and  economic  competitiveness.  The  mission  of  the  Center  is  to 
strive  for  a  productive  education  and  research  environment,  with  innovative  edu¬ 
cational  teeiching  techniques,  effective  use  of  state-of-the-art  technology,  seminars, 
visitors,  conferences,  and  publications  as  its  visible  signs  of  success.  It  was  pro¬ 
posed  that  the  CCS  become  the  primary  location  of  activities  supported  by  this 
grant.  Funding  from  this  grant  was  instrumental  in  developing  an  infrastructure 
for  research,  where  there  is 

•  release  time  for  research  faculty; 

•  summer  research  activities; 

•  up-to-date  computing  facilities; 

•  visiting  scholars’; 

•  support  for  imdergraduate  students. 
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2.  Scientific  Program 

2.1.  Optimization  and  Control  Strategies  in  Fiber  Manufacturing 

The  primary  objective  of  this  research  is  to  merge  novel  methods  for  control 
of  polymer  extrusions  with  focused  fiber  manufacturing  issues.  The  secondary 
objective  is  to  develop  finite  element  edgorithms  that  solve  specific  model  equa¬ 
tions.  When  combined  with  experiments  these  algorithms  yield  fundamental  re¬ 
lationships  between  material  properties  of  industrial  fluids  (agricultural  sprays  or 
polymeric  hquids)  and  measurable  flow  length  scales  found  in  laboratory  experi¬ 
ments.  Each  of  these  objectives  involves  the  development  and  application  of  new 
mathematical  and  numerical  methods  for  problems  defined  by  the  industrial  and 
mzuiufactiuing  community. 

In  ploymer  extrusion  processes,  one  is  often  interested  in  maintaining  a  quasi¬ 
uniform  temperature,  so  as  to  reduce  material  (extrusion  product)  inhomogeneity, 
throughout  the  extrudate  cooling  process.  As  a  preliminary  step  we  have  studied 
the  following  simphfied  case:  seek  a  uniform  temperature  distribution  at  the  exit 
\mder  a  steady  state  condition.  To  achieve  a  uniform  temperatm^  distribution  at 
the  exit  we  adjusted  the  heat  flux  on  the  smface  of  the  pipe  near  the  exit. 

The  extrudate  in  question  is  assumed  to  be  a  viscoelastic  fluid  of  the  Oldroyd 
type,  e.g.,  polymer  melts  with  a  fast  relaxation  mode.  Here  u  denotes  the  velocity 
field,  p  the  pressure  field,  T  the  temperature  field,  and  r  the  pmely  elastic  pairt 
of  the  extra  stress.  Also,  D  =|(Vu  -  Vu^),  W  =  |(Vu  -1-  Vu^)  and  DaT  = 
(u  •  V)r  +  rW  -  Wr  -  a(DT  +  tD),  where  -1  <  a  <  1.  Next,  let  Re,  We 
and  (aj  be  the  Reynolds  number,  Weissenberg  number  and  retardation  parameter, 
respectively.  With  these  definitions  the  governing  equations  for  an  Oldroyd  type 
fluid,  in  dimensionless  form,  is  given  by  the  Navier-Stokes  equations 

Re(u- V)u-|- Vp  =  (1 -a;)  Au-t-dwr-hf,  in  Q,,  (2.1) 

the  incompressibility  constraint 

div  u  =  0  in  fl,  (2-2) 

the  constitutive  equation  (Oldroyd  Model) 

T  +  We  DaT  =  2u;D  in  Q,  (2.3) 

the  boundary  condition  (for  this  prehminary  study) 

u  =  h  on  r,  (2.4) 
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the  energy  equation 


» 


— «;V7’+ (u  •  V)T  =  Q +  2/i  (Vu  + Vu^)  :  (Vu  + Vu  in  f2,  (2.5) 

and  with  additional  boimdary  conditions 

0  on  Fd,  (2.6) 

Hjv  oraFjvUro,  (2.7) 

g  on  Tc-  (2.8) 

In  this  setting  the  data  functions  f ,  Q,  Hjv,  and  h  eire  assumed  known;  the  control 
g  is  to  be  determined  so  that  hot  spots  are  avoided.  The  constants  k  and  fi  depend 
on  the  thermal  conductivity  coefficient,  density,  specific  heat  at  constant  volume, 
and  viscosity  coefficient  of  the  fiuid.  We  assumed  that  buoyancy  effects  can  be 
neglected,  and  thiis  the  temperature  variable  does  not  appear  in  (2.1). 

Two  meains  of  obtaining  a  uniform  temperature  distribution  have  been  con¬ 
sidered.  The  first  involves  making  the  gradient  of  the  temperatiu^  along  the 
boundary  To  small.  Thus,  for  example,  given  a  velocity  field  u,  we  would  seek  a 
temperature  field  T  and  control  field  g  such  that  the  functional 

=  |v,r|'  I 

is  minimized,  of  course,  subject  to  the  constraints  imposed  by  the  fiow  equations 
(2.1)  —  (2.5).  Here,  Vs  denotes  the  surface  gradient  operator,  e.g.,  in  3?^,  the 
tangential  derivative  operator  d/dr.  The  non-negative  parameters  a  and  6  can  be 
used  to  change  the  relative  importance  of  the  two  terms  appearing  in  the  definition 
of  M .  These  constants  may  also  be  used  as  penalty  parameters.  The  appearance 
of  the  control  g  in  the  definition  of  M  is  necessary  because  we  are  not  imposing 
a  priori  limits  on  the  size  of  this  control.  The  minimization  of  (2.9)  results  in  a 
quasi-uniform  temperature  distribution  along  the  boundary  segment  Fq.  This  is 
a  result  of  the  fact  that  surface  derivatives  of  the  temperature  are  forced  to  be 
small.  Another  means  of  eichieving  the  desired  result  is  to  try  to  directly  force 
the  temperature  field  itself  to  be  quasi-uniform.  In  this  setting,  given  a  velocity 
field  u,  we  would  seek  a  temperature  field  T  and  a  control  field  g  such  that  the 
functional 

3  (T, ^  -  ^■<1"  + y  /  1*1' 


dn 

dn 
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is  minimized  subject  to  (2.1)  —  (2.8),  where  is  some  desired  temjjerature  dis¬ 
tribution,  e.g.,  a  quantity  close  to  the  average  temperature  along  Fc  for  the  un¬ 
controlled  system.  Once  again  the  non-negative  parameters  7  and  6  can  be  used 
to  change  the  relative  importance  of  the  two  terms  appearing  in  the  definition  of 
S  as  well  as  to  act  a  penalty  parameters. 


2.1.1.  The  Optimization  Problem,  Existence  of  Solutions,  An  Optimal- 
^  ity  System 

Here  we  give  a  precise  statement  of  the  optimization  problem  and  results  obtained. 
We  first  note  that  (2.1)-(2.4)  uncouples  from  (2.5)-(2.8).  We  may  solve  for  (u,pt) 
from  (2.1)-(2.4)  once  and  for  edl  and  then  substitute  them  into  (2.5)-(2.8).  Thus, 
the  o^y  state  variable  is  T,  i.e.,  the  temperature  field,  and  the  only  boundary 
•  control  variable  is  g.  The  state  eind  control  variables  are  constrained  to  satisfy 

the  system  (2.5)- (2.8),  which  are  recasted  into  the  following  weak  form:  find 
T  G  =  {S  €  (fi) :  -S'  =  0  on  T^}  such  that 


«a(r,-S)-Hc(u,T,.S)  =  (Q,5')-h«(g,5')r^  +  «(H^5)r^  ^SeHl,{n),  (2.11) 
where 

Q  =  Q  -H  2/i  (Vu  -h  ( Vu)^ j  :  [Vu  +  ( Vu)^  j  . 

For  each  possible  control  function  g,  we  have  established  the  existence  of  a 
unique  corresponding  state  function  T. 


Lemma  2.1.  For  every  g  €  (Fc, ),  there  exists  a  unique  T  G  H})  (fl)  such  that 
(2.11)  is  satisfied.  Moreover, 


lini.  +  rilo,r„  <  C  (ll8llo,rc  +  IIQIIo  +  l|H-^llo,r»)  ■ 

The  admissibility  set  is  defined  by 


Uad  =  {(^,  g)  e  Hi  (n)  X  (Fc) :  S  (T,  g)  <  00,  (2.11)  is  satisfied}  . 


Then,  ^T,  G  Uad  is  call  an  optimal  solution  if  there  exists  s  >  0  such  that 


s(f',^)  <S(r,g)  V(T,g)G£/ad 
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satisfying 


T-T 


+  lie  -  g|lo,rc  ^ 


Making  use  of  the  Lemma,  we  obtained  existence  and  uniqueness  of  optimal 
solutions.  We  also  obtained  the  optimality  condition 


9  =  - 


(2.12) 


where  ^  is  the  solution  of  the  adjoint  state  equation 

Ka{R,^)  +  c  (u, i?, i  (i?, T - Td)^^  =0  (2.13) 

Elimination  g  from  (2.11)  and  combining  with  (2.13),  we  obtain  the  optimality 
system 

Ka{T,S)  +  c (u, r,  S)  +  ^(^, S)ra  =  (Q,  S)  +  k  (H^,  5)r^  V  5  G  H^n), 

(2.14) 

and 

/ca(i2,#)  +  c(u,i?,#)-i(i?,r-rd)rc=0  V  ReHliQ).  (2.15) 

o 

Thus,  the  optimal  state,  i.e.,  the  temperature  distribution  T,  can  be  found  by 
solving  the  coupled  system  (2. 14)- (2. 15),  which  also  provides  the  optimal  co-state 
The  optimal  control  g  can  then  be  deduced  from  (2.12). 


2.1.2.  Finite  Element  Approximation,  Error  Estimates,  and  Iterative 
Methods 

Next,  we  constructed  a  finite  element  algorithm  for  determining  approximations  of 
the  optimality  system  (2.14)  —  (2.15)  as  follows:  seek  G  and  G  such 
that 

Ko(r‘,S'‘)+c(u,r‘,S'‘)+^($‘,S‘)r,,  =  (Q,S'‘)+K(H„,S"‘)r„  V  S*  €  W"* 

(2.16) 

and 

Ka{R!^,^^)+c{u,R^,^^)-^{R\T^-Td)rc=0  V  R^eW^.  (2.17) 

Although  the  optimality  system  is  Unear,  the  coupling  of  u  and  ^  makes  the 
derivation  of  the  error  estimates  nontrivial.  To  obtain  our  error  estimates  we 
applied  the  Brezzi-Rappaz-Raviart  theory. 


Theorem  2.2.  Let  (T,^)  and  be  the  solutions  of  (2. 14)- (2. 15),  and 

(2.16)-(2.17),  respectively.  Assume  that  T,  ^  D  H}){Q)  for  some 

1  <m<k]  also  that  the  approximation  property 

inf  \\S  -  5'*|L  <  ChT  ||5|L .  j  V5  €  Hh{n)  andO<m<k  (2.18) 
holds.  Then, 

||T  -  T'^IIj+II#  -  <  Cmax  |i,  l}  (lIQIL-i  +  ll-ff^llr^,,m-i/2  +  ll'^<illr^„m-i/2) 

(2.19) 

where  C  is  independent  of  h,  6, 7,  T,  and 


An  iterative  algorithm  for  solving  (2.14)-(2.15)  is  defined  as  follows: 
choose 

for  n  =  1,2, ...,  solve,  for 

(T<">,  S)+c  (u,  Tl"),  S)  =  S)rc+(<3. «)+«  (H’jv,S)r»  V  S  e  HUil)-, 

then  solve,  for 


K  a  (R,  #(”))  +  c  (u,  R,  «^W)  =  i(i?,  -  Td)Tc  V  H  €  HUQ). 

The  convergence  of  this  algorithm  was  obtained  by  observing  that  it  is  effec¬ 
tively  a  gradient  method  for  the  following  minimization  problem:  find  g  G  L^{rc) 
such  that  /c(g)  =  3(r(g),g)  is  minimized  where  T(g)  G  H}){d)  is  defined  as  the 
solution  of  (2.11). 


2.1.3.  Computational  Techniques 

Computational  results  have  been  obtained  using  the  most  powerful  computational 
technique  developed  thus  far,  and  making  effective  use  of  high-performance  com¬ 
puter  technology.  In  particular,  computational  results  for  a  model  with  ein  01- 
droyd  type  fluid  have  been  obtsdned,  demonstrating  the  effectiveness  of  the  theory 
and  code  as  well  as  their  potential  applicability  to  problems  of  interest  to  the  Air 
Force  Office  of  Scientific  Research. 

We  end  this  section  by  presenting  computational  examples  considered.  Let 
be  the  unit  square  (0, 1)  x  (0, 1).  Let  F  =  rcUr'i3Ur'jvUr'o  be  as  shown  in 
figure  1. 


Figure  1.  Computational  domain 


The  finite  element  spaces  are  chosen  to  be  piecewise  linear  elements  on  a 
triangular  mesh.  All  the  numerical  results  make  use  of  the  following  parameters 
and  data: 


parameters:  k  =  1/0.73;  Re  =  1;  u  =  1/2 


ffr 

boundary  data:  T  =  1  on  Fx?;  —  =  0  onF^vj^Fo; 


dn 


heat  source 


Q  =  — Y  (  97r^cos(37rar)cos^(7ry)  —  47r^sin^(-7rx)cos(27ry)  j 


IStt 

+(1  -  y^)-^  sin(37ra:)  cos^(7ry); 


volcity  profile:  u  =  (l  — y^,0), 


elastic  extra  stress  profile: 


T  = 


/O  -y 
-y  -4y^ 


pressure  profile:  p  =  A{\  —  y^)  —  x. 


In  the  functional  (2.10)  we  choose 


For  the  data  given  above,  the  exact  solution  of  the  uncontrolled  problem,  i.e.. 


for 


dn 


on  Fc, 


is  given  by  T  =  5sin^(|7ra:)  cos^(7r^)  +  1. 

We  compare  the  temperature  distribution  in  the  uncontrolled  case  with  the 
optimal  temperature  distribution  in  the  controlled  case  for  which 


dn 


9  on  Tc, 


where  g  is  the  control  such  that  the  functional  (2.10)  is  minimized.  Approxima¬ 
tions  to  the  optimal  state  and  co-state  are  computed  form  (2.14)-(2.15);  the  ap¬ 
proximate  optimal  control  is  then  obtained  from  (2.12),  i.e.,  g'^  =  — (1/5)4*  |r^ . 
AU  of  the  computational  results  shown  below  were  obtained  with  the  use  of  mesh 
size  h  =  ^.  Of  course,  calculations  with  varying  mesh  sizes  were  performed.  Since 
these  merely  verified  the  errors  estimates,  we  do  not  report  on  them  here. 

Specifically,  figures  2-4  deal  with  the  following  cases: 

2.  Exact,  uncontrolled  temperature  contours  £ind;  Optimal  temperature; 


Figure  2.  Temperature  contours:  uncontroled  and  controled. 
Optimal  boundary  control  on  Tc- 
(Fc  is  the  top  boundary  segment.) 

(Fo  is  the  right  boundary  segment.) 
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3.  Exact,  uncontrolled  temperature  smfaces  and;  Optimal  temperature; 
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# 


# 


2.2.  Simulation  of  Maguetoelastic  Interactions 

Terfenol-D,  which  is  Tbj.  Dyi_x  Fe^,  x  ^  0.3,  is  a  magnetoelastic  material  ex¬ 
hibiting  a  complicated  microstructure  that  responds  easily  and  drcunatically  to 
an  applied  magnetic  field.  It  has  the  largest  known  magnetoelastic  coupling,  with 
huge  magetostriction  and  small  magnetomechanical  losses,  giving  it  extraordi¬ 
nary  promise  for  use  as  an  actuator  and  sensor  material.  Its  microstructure  or 
fine  scale  morphology  has  a  clear  role  in  its  macroscopic  behavior.  Two  issues 
are  paramount.  First,  the  presence  of  several  stable  variants  at  a  given  temper¬ 
ature  refiects  a  complicated  potential  weU  structure  for  the  free  energy  of  the 
materiad.  Second,  the  presence  of  spatially  oscillatory  behavior  at  the  small  scale 
suggests  competition  between  the  free  energy  of  the  material  and  loading  or  other 
environmental  effects.  Both  of  these  features  represent  highly  nonlinear  processes. 

For  large  specimens,  such  as  those  encoimtered  in  usual  actuator  applications, 
the  presence  of  complicated  domain  structures  has  hindered  the  usefulness  of 
theories  which  seek  to  incorporate  elastic  eflfects.  A  theory  of  magnetostriction 
which  is  micromagnetic  in  natiue  has  been  developed  by  consulteints  Richeird 
James  and  David  Kinderlehrer.  The  focus  is  on  the  precise  specification  of  the 
potential  wells  of  the  anisotropy  energy.  There  are  many  parallels  between  shape- 
memory  and  giant  magnetostrictive  behavior.  A  key  difference  is  that  shape 
memory  materials  undergo  a  first  order  transition  resulting  in  2-15%  streiin  while 
Terfenol  suffers  a  second  order  transition  with  only  .25%  strain. 

By  advancing  and  implementing  analytical  methods  it  has  been  shown  that  the 
theory  of  magnetostriction  successfully  predicts  the  equilibrium  microstructures 
observed  in  Terfenol-D.  The  theory  furthermore  suggests  a  mechsinism  for  the 
giant  magnetostriction  related  to  exchange  of  stability  of  the  variants  within  a 
laminate  or  among  laminate  systems  on  which  a  magnetic  field  has  been  applied. 

To  understand  this  phenomenon  more  completely,  in  particular  away  from 
equilibrium,  the  interdisciplinary  group  has  established  a  broad  program  com¬ 
posed  of  experiment,  theory,  and  simulation. 

2.2.1.  Simulation  Program 

The  energy  of  the  magnetoelastic  systems  is  described  in  terms  of  a  free  energy 
of  the  form 

£'(y,m,H)  =  j  {lF’(Vy(x),m(y(a:),r)  —  H  •m(y))det  Vy(a;)}  dx  +  f  iVttf  dx 
Jn  Jm 

(2.20) 
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|det  Vy(a?)m|  =  1  in  and  Au  =  div  (mXj,(o))  in  9?^, 

where  y{x)  is  the  elastic  deformation  ,  C  is  a  reference  domain,  m(y)  the 
magnetization  ,  and  H  the  apphed  magnetic  field.  The  second  equation  in  (2.20) 
represents  the  magnetostatic  equations  from  Maxwell’s  Theory.  Since  we  are 
dealing  with  a  large  body  approximation,  exchange  energy  is  ignored.  The  stored 
energy  density  W (A,  m,  x)  is  highly  nonconvex  and  gives  rise  to  a  functional  which 
is  not  lower  semicontinuous. 

The  simulation  program  started  with  the  simulation  of  hysteretic  behavior  in 
two  dimensions  for 

£'(H,in)  =  f  (0(m,a;)  —H  *111)010;  +  ^  f  |Vit|^  dx  (2.21) 
Ja  ^  Jm 

Au  =  div  mxn  in  |in|  =  1  in  (2.22) 

where  an  applied  field  is  introduced  in  the  system  and  cycled  through  a  range  of 
values. 

Using  methods  of  nonlinecu  analysis  computational  results  and  diagnostics 
have  been  explored.  Since  the  energy  picture  is  mesh  independent,  computing  on 
a  fairly  coarse  grid  sufiices  to  estabhsh  its  character.  Kinderlehrer  £ind  Ma  deter¬ 
mined  analyticedly  the  hysteresis  loop.  For  a  magnetic  system,  this  analysis  rests 
on  the  introduction  of  a  shadow  energy  for  a  simplified  version  of  the  system. 
It  is  believed  that  this  simplified  version  may  suggest  possible  dispersive  inter¬ 
actions  which  may  be  attributed  to  a  shape-memory  or  pseudoelastic  body.  It 
was  observed  that  the  simulation  displays  a  modified  Stoner- Wohlfarth  scenario. 
This  means  that  the  system  evolves  through  stable/metastable/unstable/stable 
regimes.  By  modeling  the  nature  of  the  instability  one  is  able  to  determine  the 
critical  field  event. 

These  collaborators  than  modified  the  simulation  to  include  magnetoelastic 
effects,  in  the  context  of  linear  elasticity.  The  linear  elastic  magnetostrictive 
energy  has  the  form 


«/>(£,  a) 
e 


^  (Vu  +  Vu^) ,  m  =  Ma,  jal  =  1, 


where 


14 


m 


^elastici.^^ 


0an(a) 


2  ^  +  Ci2(s^1i£^22  +  ^22^^33  +  £332^11 )  +  2^44(2^12  +  ^23  +  ^13)5 

6otr  £  +  bi  Suaj  +  62(s^uai«2  +  £230:20:3  +  £13010:3), 

K-£at. 


A  two  dimensional  model  was  obtained  by  projecting  onto  the  (0  —  11)  plane, 
making  the  assumption  that  (peiastic  is  isotropic,  and  changing  coordinates  so  that 
(—211)  is  xi  and  (111)  is  X2.  In  this  setting 


0+(£,  a)  =  <j)aastic{^)  +  o)  +  0o„(o) 


with 


<l>eiasuM  =  /4  2^£■,•  +  ^A(tre)^  /1>0,  A>0, 

o)  =  b{£n  -  7ei2)((oi)^  -  7Q:ia2)  -  b'  Y^eijaiaj, 

‘/'^(o)  =  «((oi)^  -  7Q;l02)^  «  >  0. 

Here  7  =  2v^,  6  <  0,  6*  >  0,  with  b  approximately  equal  to  —62,  and  b'  is  roughly 
61  —  i>2*  To  calculate  the  magnetostriction  in  the  original  (—211)  direction,  now 
the  (1,0)  direction,  one  must  first  compute  the  average  of  e  and  then  apply  the 
formula.  Thus,  one  has 


=^11  • 

Unfortunately,  the  above  simulations  do  not  give  satisfying  results,  i.e.,  the 
energies  described  above  do  a  poor  job  of  resolving  the  microstructure.  The 
challenge  facing  the  FAMU  group  was  to  develop  a  better  simulation. 

2.2.2.  Computational  Techniques  and  Code  Development 

The  global  infinmim  of  the  energies  described  in  section  2  are  not  attained.  How¬ 
ever,  it  is  assumed  that  the  underlying  physical  microstructure  can  be  represented 
by  suitable  minimizing  sequences.  In  order  to  develop  a  computer  simulation 
which  can  provide  us  with  an  efiicient  and  effective  selection  principle  for  such 


minimizing  sequences,  our  strategy  is  to  employ  the  most  powerful  computational 
techniques  developed  thus  far,  develop  new  ideas  and  make  effective  use  of  high- 
performance  computer  technology.  Thus,  we  have  simultaneously  been  working 
on  two  different  approaches:  (i)  a  global  minimization  method  based  on  the  idea 
of  simulated  annealing;  and  (ii)  a  study  of  the  minuxuzation  of  a  relaxed  problem 
where  generalized  solutions  in  terms  of  Yoimg  measures  for  the  original  noncon- 
vex  problem  are  approximated.  Both  of  these  approaches  must  somehow  address 
the  following  challenging  problem:  loss  of  uniqueness  of  global  minimum,  the 
emergence  of  abundant  local  minimum  or  even  saddle  points. 

A  state-of-the-art  computer  model,  developed  by  Carnegie  Mellon  researchers, 
that  can  simulate  unconstrained  minimization  of  the  angular  variable  and  the 
displacement  for  lineeir  magnetostriction  has  been  developed.  The  hysteresis  dia¬ 
gram  is  computed  usmg  a  continuation  of  solutions  with  respect  to  increasing  and 
decreasing  the  applied  magnetic  field  along  the  xi-axis.  Work  is  ciurently  under 
way  to  extend  this  code  to  nonlinear  magnetostriction. 

Eklucational  activities  in  this  area  of  research  included  the  training  of  a  leu-ge 
number  of  minority  undergraduate  and  graduate  students.  We  have  also  offered  a 
special  lecture  series  on  modeling  and  simulation  of  problems  cirisiag  in  the  study 
of  materials  science.  Also,  special  discussion  groups  were  held.  These  activities 
brought  some  of  the  researchers  and  their  research  into  the  classroom. 

2.3.  Controllability  for  Parabolic  Equations 

In  many  problems  of  interest  to  materials  science  and  engineering  there  arises  the 
need  of  controllability  for  parabolic  equations.  The  purpose  of  this  research  was  to 
investigate  exact  controllability  problems  for  systems  modeled  by  linear  parabolic 
differential  equations.  Our  aim  of  investigation  has  been  to  examine  the  limit 
behavior  of  the  correspondiug  optimal  control  problem.  Using  this  technique  we 
obtained  a  new  condition  on  the  exact  controllability  of  the  system. 

We  now  give  a  statement  of  the  problem.  Let  be  a  bounded  domain  in 
For  fixed  T,  let  Q  =  [0,7]  x  and  E  =  (0,T)  x  dD..  Let  A  be  a  second  order 
elliptic  differential  operator  defined  by 

A{y)  =  -  XI  ^  ^  ^  y- 

ij=l  '  3  / 

We  assume  that  c>  0  on  and  that  the  matrix  (ajj  (x))is  symmetric  and  positive 
definite. 


With  these  definitions,  we  have  studied  the  following  initial  boundary  value 
problem 


■^  +  A(y)  =  0  inQ, 
^  =  u{t,x) 


on  S, 


(  y(x,0)  =  ^(x)  in 


(2.24) 


where 


= 


with  Hi  being  the  component  of  the  outward  unit  normal  to  dfl.  We  say  that 
the  system  is  exactly  controllable  for  a  given  state  y  if  there  exists  u  G  (S)  such 
that  the  solution  y  of  problem  (2.23  )  also  satisfies 

y(Tr)  =  y{-). 

Next,  we  define  the  cost  functional  Je  by 

Js{u)  =  ^  y (y(^>x,«)  -  y(x))^  dx+l  J u^{t,x.)  dxdt, 


where  y  =  y{t,x-,u)  satisfies  (2.24).  Ue{t,x)  is  said  to  be  an  optimal  control  if 
it  satisfies 

Jei'^e)  =  inf  Je{u). 

ueL^Q) 

Our  weak  formulation  of  problem  (2.24)  takes  the  form 


(^,</*) +a(y,<^)  = 

y(0,x)  =  yo{x), 


for  0  €  (fi) , 


where  the  bilinear  form  a  :  (fl)  x  (fl)  — ^  R  is  given  by 

=  y  I  ^  c(x)(/>^|  dx. 


(2.25) 
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where  A,  B,  U,  and  V  are  the  appropriately  chosen  admissible  sets. 

Remark  1.  State  equations  (2.26)  and  (2.27)  are  given  in  an  abstract  form 
that  may  represent  partial  diflferential  equations  or  more  general  functional  relar 
tions.  However,  we  shall  always  take  them  to  represent  partial  differentied  equa^ 
tions.  These  equations  are  also  written  in  strong  form.  In  order  to  place  them 
in  an  abstract  weak  formulation,  let  W  denote  a  normed  space  of  test  functions 
w(x)  =  (wi,W2),  where 

wi  e  (Lg(^l)Y,  W2  €  (L,(Q))"'  and  ^  J  =  1 
with  1  <  p  <  oo,  and  1  <  g  <  oo.  We  define  the  operator  h(w,u,u) :  W  — >  3?  by 

h{w)=  /  fo(w,u{x),M(x,u))wi  dx  +  /  go(x,{K(^,v)){x),v{x))w2  d/x 
Jn  Jt 

/ 

where  foWi  and  go'W2  denote  inner  products  in  3?’’  and  3?^  ,  respectively.  Instead 
of  state  equations  (2.26)  and  (2.27)  we  may  consider  a  single  state  equation  in  the 
weak  form 


/(w)  =  h{w), 

where  the  operator  I  repleices  both  operators  L  and  I. 

Thus,  when  considering  control  strategies,  we  are  interested  in  the  problem  of 
minimizing  functionals  of  the  form 

^{w,u,v)=  [  f{w,u{x),M{x,u))dx+  I  g{x,{K{^,v)){x),v{x))  dfi  (2.30) 
Jn  Jt 

with  state  equation  (in  weak  form) 

Z(w)  =  h{w)  for  aU  w  G  W,  (2-31) 

and  constraints 

M{XyU)  G  A,  u{x)^U  a.e.  in  fZ,  (2.32) 

K{x,u)  G  B,  u{x)gV  p  —  a.e.  on  (2.33) 


Remark  2.  Of  particuleir  interest  is  the  case  where  the  element  w  of  an 
admissible  triple  (w,u,v)  uniquely  determines  the  controls  u  and  v. 

Remark  3.  In  practice  one  may  define  the  control  set  U  by 

U={uGL°^[0,T\:  |HL<M}  (2.34) 

for  M  >  0. 

Control  Strategy  1.  A  special  case  of  the  above  abstract  optimal  control 
problems  is:  How  do  we  choose  a  control  function  g  ,  so  that  a  crystal  interface 
will  approximate  a  prescribed  desired  profile? 

Let  f2  be  a  region,  occupied  by  the  crystal.  Then,  a  mathematical  description 
of  this  control  problem  proceeds  as  follows.  If  C{x,t)  is  a  positive  continuous 
function,  and  ^(t)  is  a  continuous  positive  function  for  0  <  t  <  T.  Then,  a 
possible  cost  functional  may  be  one  of  the  form 

9(u,  g)  =  J(u)  +  ~  f  |v.9|"  dr,  (2.35) 

in  which  some  desired  optimization  objective  J(u)  is  considered  in  conjimction 
with  the  cost  of  applying  the  required  control  g,  i.e.,  the  norm  of  g.  Here  Vs  is 
the  surface  gradient  operator  on  the  interface  T  and  5  is  a  fixed,  sufficiently  small 
parameter.  The  desired  optimization  objective  J(u)  that  one  could  consider  is 
the  tracking  functional  for  the  interface.  This  functional,  with  its  domain  being 
the  interface,  has  the  form 

J(u)  =  A  f  ilr{t)[H{r,t)  —  C{r,t)Ydil.  (2.36) 

Where  A  is  a  suitable  positive  constant,  which  reffects  the  different  needs  of  each 
control  strategy. 

We  wish  to  minimize  9(u)  subject  to  some  appropriate  state  equations  and 
boundary  conditions. 

Control  Strategy  2. 

In  free-boundary  problems  one  must  use  adaptive  meshing  that  tracks  the 
interface  motion.  The  difficulty  is  that  these  schemes  may  not  be  able  to  han¬ 
dle  topological  changes  without  user  intervention.  This  in  turn  makes  real-time 
control  inconvenient  during  the  simulation.  Berg,  Yezzi,  and  Tannenbaum  have 
explored  aai  alternative  approach  [10].  In  their  approach  one  begins  by  writing 


down  the  equations  of  motion  of  the  evolving  interface  directly.  In  two  dimensions, 
these  are  the  equations  of  curve  evolution.  In  this  setting  the  moving  interface 
is  described  by  a  family  of  parameterized  curves,  C  :  [0,1]  x  [0,t/]  The 

evolution  of  the  cmve  describing  the  interface  is  governed  by 

^  =  a{s,t)T  +  l3{s,t)N  (2.37) 

where  s  parametrizes  the  curve,  N  is  the  normal  vector,  T  is  the  tangent  vector, 
and  a  and  P  are  velocity  functions.  Ciurently  we  are  only  considering  the  shape 
of  the  interface,  thus  we  may  take  a  =  0.  Note  that  changing  a  only  changes  the 
curve’s  pareimetrization,  and  not  its  shape.  In  addition,  if  we  assume  that  the 
motion  of  the  curve  is  determined  solely  by  the  local  geometry  of  the  curve,  then 


dt 


=  P{k)N 


(2.38) 


where  «(s,t)  is  the  curvature. 

There  are  several  difficidties  that  one  may  encounter  when  seeking  numerical 
solutions  of  curve  evolution  problems.  For  example,  if  curvature  terms  are  absent 
or  small,  the  solution  will  typically  develop  comers  or  shocks,  even  if  the  initial 
data  is  smooth.  Osher  and  Sethian  [11-13],  and  their  coworkers,  have  developed  a 
series  of  algorithms  that  successfully  address  many  of  the  problems  encountered 
in  curve  evolution  problems.  These  algorithms  employ  a  level  set  representation 
of  the  interface,  and  are  the  foundation  of  the  numerical  approach  developed  by 
Berg,  Yessi,  and  Tannenbaum  [14]. 

The  idea  of  applying  the  theory  of  curve  evolution  for  modeling  crystal  evolu¬ 
tion  has  been  considered  by  Srinivasan,  et  al.  [15] .  Similar  real-time  interface  esti¬ 
mation  and  control  techniques  have  been  developed  for  reactive  etching  processes 
by  a  number  of  authors.  For  example,  Shaqfeh  and  Jmgensen  [16],  and  Singh,  et 
al.  [17]  have  done  pioneering  work  in  this  area.  Adalsteinsson  and  Sethian  build 
and  greatly  extended  these  techniques  to  a  novel  approach  for  modehng  etching, 
deposition,  and  lithography  fabrication  processes  [  13,18]. 

It  is  our  intension  to  explore  the  possible  use  of  these  algorithms  for  tracking 
the  evolution  of  the  crysteil  interfaces.  The  focus  of  our  study  has  been  the 
development  of  techniques  that  directly  control  the  evolving  features  of  the  crystal 
interfaces.  This  is  a  challenging  task,  however  the  potential  benefits  are  large.  The 
level  set  techniques  show  potential  for  real-time  control. 
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2.4.  Controllability  for  Parabolic  Equations 

In  many  problems  of  interest  to  materials  science  and  engineering  there  arises  the 
need  of  controllability  for  parabolic  equations.  The  purpose  of  this  research  was  to 
investigate  exact  controllability  problems  for  systems  modeled  by  linear  parabolic 
differential  equations.  Our  aim  of  investigation  has  been  to  examine  the  limit 
behavior  of  the  corresponding  optimal  control  problem.  Using  this  technique  we 
obtained  a  new  condition  on  the  exact  controllability  of  the  system. 

We  now  give  a  statement  of  the  problem.  Let  be  a  bounded  domain  in 
For  fixed  T,  let  Q  =  [0, 7^  x  and  S  =  (0,  T)  x  Let  A  be  a  second  order 
elliptic  differential  operator  defined  by 

My) = -  ^  ^  ^^)  y- 

ij=l  '  3  / 

We  assume  that  c  >  0  on  and  that  the  matrix  (ajj  (x))is  symmetric  and  positive 

definite. 

With  these  definitions,  we  have  studied  the  following  initial  boimdary  value 
problem 

%^A(y)  = 

^  _ 

du 

(  2/(x,0)  = 

where 

a  A  d 

with  rii  being  the  component  of  the  outward  unit  normal  to  dO,.  We  say  that 
the  system  is  eocactly  controllable  for  a  given  state  y  if  there  exists  u  ^  (S)  such 

that  the  solution  y  of  problem  (2.39  )  zilso  satisfies 

y(T,-)  =  v(-)- 


0  in  Q, 

u(t,  x)  on  S, 

^(x)  in  (2.40) 


Next,  we  define  the  cost  functional  Jg  by 


Je{u) 


\  J {yiT,x,u)  -  y{-x.)f  dx+l  J w^(i,x)  dxdt, 
n  s 


where  y  =  y{t,x\u)  satisfies  (2.40).  Us{t,x)  is  said  to  be  an  optimal  control  if 
it  satisfies 

Our  weak  formulation  of  problem  (2.40)  takes  the  form 

for  e  fl-' (!J) . 

y(0,x)  =  yo(x),  (2.41) 

where  the  bilinear  form  a  :  (fi)  x  (fi)  — >  R  is  given  by 


a(0,'0) 


aij(x) 


d4>  dip 

dxi  dxj 


+  c{x)4np 


) 


dx. 


Our  main  results  of  this  research  has  been  obtained  by  using  an  operator 
representation  of  the  terminal  state.  That  is,  we  define  an  operator  R  :  L^{Q) 
L^(f2)  by 

(s  +  R){y(T)-y)=e(y-E{T)v) 
where  E{t)  is  the  solution  operator  for  the  problem 

on  Q,  (2.42) 

on  E,  (2.43) 

i«(0,x)  =  2:(x)  in  Cl,  (2.44) 

i.e.  w{t,x)  =  E{t)z{x). 

The  following  properties  of  the  operator  R  have  been  established:  R  is  a 
compact  operator  in  L^(fl);  R  is  symmetric  and  semi-definite;  if  S  is  analytic, 
then  ker  R  =  Q. 


dw 

dv 


0 

0 
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Note  that  system  (2.42)  -  (2.44)  is  said  to  be  approximately  controllable  if,  for 
any  given  y  and  5  >  0,  there  exists  a  control  u{t,  x)  and  a  function  y{t,  x)  such 
that  (y(t,x),u(t,x))  is  a  solution  of  (2.40)  aad 

||y(T,x)-y(x)||  <8. 

Approximate  controllability  is  weU  known  for  linear  parabolic  differential  equa¬ 
tions.  However,  the  proof  that  we  were  able  to  obtain  is  different  from  standmd 
proofs,  in  that  it  is  constructive. 

Finally,  using  properties  of  R  we  obtained  the  following  result  on  exact  con¬ 
trollability  of  system  (2.40). 

Theorem  2.3.  Let  he  the  eigensystem  of  R.  If 

y-E{T)yoeD{Rr^), 
then  system  (2.40)  is  exactly  controllable. 
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Temperature  Control  In  Polymer  Extrusion  Processes! 


Max  D.  Gunzburger 
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James  C.  Turner 


Abstract.  An  optimization  problem  is  formulated  motivated  by  the  desire  to  obtain  uniform  extrudate  tcmpoiature  at 
the  die  exit  in  a  polymer  extrusion  process.  Control  is  effected  by  adjustments  to  the  heat  flux  along  the  surface  of  the  pipe. 
An  optiniality  system  of  partial  differential  equations  is  derived  from  wliich  optimal  controls  and  states  may  be  deteriiiincd. 
Then,  finite  clement  discretizations  of  the  optimality  system  are  defined  and  error  estimates  are  provided  along  with  an  efficient 
5oliilii>n  algorithiu  for  the  discrete  equations.  Finally,  computational  results  are  given  for  a  model  example  with  Olthiiyd  type 
fluid,  demonstrating  the  elfectiveness  of  our  theory  and  method  as  well  as  their  potential  applicability  to  industrial  pn>blems. 

1.  liiirodiictiou 

In  polymer  extrusion  processes,  one  is  often  interested  in  maintaining  a  quasi-uniform  temperature,  so 
as  to  reduce  material  (extrusion  product)  inhomogeneity,  throughout  the  extrudate  cooling  process.  In  this 
paper,  we  will  study  the  somewhat  simplified  case:  one  tries  to  obtain  a  uniform  temperature  distribution  at 
the  exit  under  steady  state  situation.  The  means  we  use  to  achieve  such  a  uniform  temperature  distribution 
at  the  exit  is  to  adjust  the  heat  flux  on  the  surface  of  the  pipe  near  the  exit. 

The  extrudate  in  question  will  be  assumed  to  be  viscoelastic  fluid  of  Oldroyd  type  such  as  polymer  melts 
with  a  fast  relaxation  mode  [10].  Let  u  denote  the  velocity  field,  p  the  pressure  field,  T  the  temperature  field 
and  T  tlie  purely  elastic  part  of  the  extra  stress.  Let  D  =  i (grad  u  +  grad  W  =  ^(grad  u  -  grad  u'^') 
and  DaT  =  (u  grad  )t+ rW  -  Wt- a(DT+ rD)  where  -1  <  a  <  1.  The  parameters  Re,  We  and  u  are  the 
Reynolds  iiiiinber,  VVeissenberg  number  and  retardation  parameter,  respectively.  The  governing  equations 
fur  the  Oldroyd  type  fluid,  in  dimensionless  form,  is  given  by  the  Navier-Stokes  equations 

/ie(u  •  grad )u  +  grad p  =  (1  —  w)Au  +  divr  + f ,  in  Q  ,  (1.1) 


the  incompressiblity  constraint 


divu  =  0  in  . 


(1.21 


the  constitutive  equation  (Oldroyd  model) 


T+WeDaT=2uT>  infi, 
and,  for  simplicity,  the  boundary  condition 


u  =  h  on  r , 


(1.3) 


(1,4) 


and  also  the  energy  equation 

-/cAT+(u-grad)T 

=  Q  +  2;i(grad  u  +  grad  u’’)  :  (grad  u  +  grad  u^)  in  fi  , 


(1.5) 


with  boundary  conditions 

T=0  onFi,,  (1-6) 

^  =  Hn  onr^UTo,  (1-7) 

on 

^=g  on  Tc.  (1-8) 

on 

The  data  functions  f ,  Q,  Hjf,  and  h  are  assumed  known;  the  control  ^7  is  to  be  determined  so  that  hot  spots 
are  avoided.  The  constants  k  and  fi  depend  on  the  thermal  conductivity  coefficient,  density,  specific  heat 
at  consteint  volume,  and  viscosity  coefficient  of  the  fluid.  See  [11]  for  details.  We  assumed  that  buoyancy 
effects  can  be  negelected,  and  thus  the  temperature  variable  does  not  appear  in  (1.1). 

Two  means  of  obtaining  a  uniform  temperature  distribution  come  to  mind.  The  first  is  to  make  the 
gradient  of  the  temperature  along  the  boundary  T,  small.  Thus,  for  example,  given  a  velocity  field  u,  we 
would  seek  a  temperature  field  T  and  a  control  field  g  such  that  the  functional 

M(T,g)  =  ^  f  [grad,  Tp  dr  +  ^  /  |<7p  dT  (1.9) 

^  Jro  ^  •'rc 

is  minimized  subject,  of  course,  to  the  constraints  imposed  by  the  flow  equations  (1.5)-(1.8).  Here,  grad, 
denotes  the  surface  gradient  operator,  e.g.,  in  R^,  the  tangential  derivative  operator  d/dr.  The  non-negative 
parameters  a  and  6  can  be  used  to  change  the  relative  importance  of  the  the  two  terms  appearing  in  the 
definition  of  A4  as  well  as  to  act  as  penalty  parameters.  The  appearance  of  the  control  g  in  the  definition  of 
is  necessary  because  we  are  not  imposing  any  a  priori  limits  on  the  size  of  this  control.  The  minimization 
of  (1.9)  results  in  a  quasi-uniform  temperature  distribution  along  the  boundary  segment  because  the 
surface  derivatives  of  the  temperature  are  forced  to  be  small.  Another  means  of  achieving  the  desired  result 
is  to  try  to  directly  force  the  temperature  field  itself  to  be  quasi-uniform.  Thus,  now,  given  a  velocity  field 
u,  we  would  seek  a  temperature  field  T  and  a  control  field  g  such  that  the  functional 

=  \T-T,\^dr  +  !^  I  lypdr  (i.io) 

is  minimized  subject  to  (1.1)-(1.8),  where  Td  is  some  desired  temperature  distribution,  e,g.,  something  close 
to  the  average  temperature  along  Tc  for  the  uncontrolled  system.  The  non-negative  parameters  7  and  6  can 
be  used  to  change  the  relative  importance  of  the  the  two  terms  appearing  in  the  definition  of  *7  as  well  as  to 
act  as  penalty  parameters.  As  will  be  demonstrated  by  numerical  examples  in  §6,  a  small  7  is  more  useful 
in  achieving  quasi-uniform  boundary  temperature  distributions,  although  it  also  reduces  the  accuracy  of  the 


approximate  solution.  We  will  examine  the  latter  issue  in  §4.  Numerical  experiments  show  that  both  (1.9) 
and  (1.10)  work  effectively  for  the  desired  objective.  We  will  focus  our  discussion  on  (1.10)  thoroughout  this 

paper. 

Under  the  realistic  assumption  that  u  -  n  =  0  on  Tc  U  Tiv,  we  may  prove  the  existence  and  uniqueness 
of  optimal  solutions  and  derive  an  optimality  system,  i.c.,  a  set  of  equations  from  which  the  optimal  control 
and  state  may  be  determined.  Also,  finite  element  methods  are  used  to  compute  an  approximate  solution 
of  the  optimality  system.  Optimal  error  estimates  are  derived  and  numerical  examples  are  presented.  We 
have  also  developed  an  iterative  algorithm  to  compute  the  approximate  solution.  The  convergence  of  our 
algorithm  is  proved  and  a  comparison  with  the  direct  method  is  made. 

Control  problems  for  the  fully  coupled  problem  as  well  as  temperature  matching  for  the  entire  extrudate 

will  be  addressed  elsewhere. 

We  close  this  section  by  introducing  some  of  the  notation  used  in  subsequent  sections.  Throughout,  C 
will  denote  a  positive  constant  whose  meaning  and  value  changes  with  context.  Also,  iy*(2)),  s  G  R ,  denotes 
the  standard  Sobolev  space  of  order  s  with  respect  to  the  set  D,  where  V  is  either  the  flow  domain  fi,  or  its 
boundary  F,  or  part  of  that  boundary.  Of  course,  H^{V)  =  V^{V),  Dual  spaces  will  be  denoted  by  (•)  .  Of 
particular  interest  will  be  the  space 

H}){^l)  =  {SeH\Q)  :  5  =  0  onFp}. 

Norms  of  functions  belonging  to  H^(Q)y  H^(T)  and  H^(Tc)  denoted  by  ||  *  i  ||  •  ||3,r  ^.nd  ||  •  Ha.rc) 
respectively. 

The  inner  product  in  L^(Q)  is  denoted  by  (•,*))  ^  (‘jOr*  in  L^{To)  by 

that  in  L^iTc)  by  (-jOrc*  Since,  in  general,  we  will  use  i^-spaces  as  pivot  spaces,  these  notation  will  also 
be  employed  to  denote  pairings  between  Sobolev  spaces  and  their  duals. 

We  will  use  the  bilinear  form 

a(r,5)=  /  gradT  gradSdn  WT,SeH^{n) 

Jfi 

and  the  trilinear  form 

c(u,r,5)  =  /(u.gradr)5<in  Vu€H*(fi)  and 

Jn 

These  forms  are  continuous  in  the  sense  that  there  exist  constants  Ca  and  Cc  >  0  such  that 

KT,5)|  <  c«||Tl|i|151h  (1.11) 

and 

|c(u,T,5)|<c,||u||i||T||i||5|li  Vu€H^(fi)  and  (1-12) 

Moreover,  we  have  the  coercivity  property 

a{T,T)  >  C<:\\T\\l  (1-13) 

for  some  constant  Ca  >  0. 

For  details  concerning  the  notation  employed  and  the  inequalities  (1.11)- (1.13),  one  may  consult,  c.p., 
[1]  and  [7]. 


2.  The  Optimization  Problem,  Existence  of  Solutions,  And  Optimality  System 

We  begin  by  giving  a  precise  statement  of  the  optimization  problem  we  consider.  We  will  assume  the 
domain  0  is  a  polygon  in  R*.  We  first  recall  that  (1.1)-(1.4)  uncouples  from  (1.5)-(1.8).  We  may  solve  for 
(u,pr)  from  (l.l)-(l-4)  once  and  for  aU  and  then  plug  them  into  (1.5)-(1.8).  Thus  the  only  state  variable  is  T, 
i.e.,  the  temperature  field,  and  the  only  boundary  control  variable  is  g.  The  state  and  control  variables  are 
constrained  to  satisfy  the  system  (1.5)-(1.8),  which  we  recast  into  the  following  weak  form:  find  T  € 
such  that 

Ka{T,S)  +  c{n,T,S)  =  {Q,S)  +  K{g,S)rc  +  K{HN,S)rN  '^SeH}}{Q),  (2.1) 

where  we  have  introduced  the  simplifying  notation 

Q  =  Q  +  2p(grad  u  +  grad  u^)  :  (grad  u  +  grad  u^) . 

Note  that  since  we  seek  T  £ 

T  =  0  on  To-  (2-2) 

Throughout,  we  will  assume  that  the  given  velocity  field  u  is  smooth  and  satisfies 

divu  =  0  in  n  and  un>0  a.e.  on  Fc  U  F^v .  (2.3) 

Under  these  assumptions,  we  have  the  useful  relation 

c(u,5,5)  =  i  / (u  11)52 dF  >0  V5  €  H^Q) , 

2  Jr 

which  can  be  derived  by  setting  T  =  5  in  the  following  integration  by  parts  formula: 

c(u,  T,  S)  =  Jjiu  ■  n)TS  dT  -  c(u,  5,  T) . 

For  each  possible  control  function  g,  there  exists  a  unique  corresponding  state  function  T. 

Lemma  2.1-  For  every  g  €  L'^{Tc),  there  exists  a  unique  T  €  such  that  (2.1)  is  satisfied 

||r||i  +  ||T|lo,rc  <  C'diffllo.rc  +  IIQIIo  +  II^Ar||o.rAr)  • 

The  admissibility  set  Uad  is  defined  by 

Uad  =  {(T,  g)  €  X  L\Tc) :  J(T,  g)  <  00,  (2.1)  is  satisfied}  .  (2.6) 

Then,  (f,g)  £  Uad  is  called  an  optimal  solution  if  there  exists  c>  0  such  that 

J{T,  9)  <  g)  V  (T,  g)  £  Uad  satisfying  ||r  -  f\\i  +  Uff  -  ff||o,rc  <  « •  (2-7) 

Based  on  the  previous  lemma,  we  can  show  the  existence  and  uniqueness  of  optimal  solutions. 
Theorem  2.2-  There  exists  a  unique  optimal  solution  {f,g}  £Uad-  ® 

Using  techniques  in  e.g.  [9]  we  may  obtain  the  optimality  condition 


(2.4) 


Moreover, 

(2.5)1 


where  $  is  the  solution  of  the  adjoint  state  equation 


Ka(«,$)  +  c(u,iz,$)-i(iz,r-rd)r^=o  \/ReHhin).  (2.9) 

Eliminating  g  from  (2.1)  and  combining  with  (2.9),  we  obtain  the  opUmalHy  system 

K a{T, S)  +  c(u,  T, S)  +  ^($, 5)rc  =  (Q.  S)  +  k  (ffjv, 5)r^  V  S  €  (2.10) 

and  . 

Ka(R,^)  +  c(u,R,<^)--(R,T-Td)r^  =0  V  R  €  RoW  ■  (2-11) 

Thus,  the  optimal  state,  t.e.,  the  temperature  distribution  T,  can  be  found  by  solving  the  coupled  system 
(2.10)-(2.11),  which  also  provides  the  optimal  co-state  The  optimal  control  g  can  then  be  deduced  from 

(2.8). 


3.  Finite  Element  Approximation,  Error  Estimates,  and  Iterative  Methods 

In  the  usual  manner,  one  may  construct  finite  element  subspaces  C  nC(fl)  parametrized 

by  a  parameter  h  that  tends  to  zero.  (In  practice,  h  is,  of  course,  related  to  a  grid  size.)  We  assume  the 
approximation  property  (see  [3]):  there  exist  an  integer  k  and  a  constant  C  such  that 


inf  ||5-5'‘||i<C'h”*||5||„.+i  VS  (n)  and  0<m<l:.  (3.1) 

A  finite  element  algorithm  for  determining  approximations  of  the  solution  of  the  optimality  system  (2.9)- 
(2.10)  is  then  defined  as  follows:  seek  T'*  €  W'  and  G  W*  such  that 


K  a{T^ ,  S'*)  -I-  c(u,T\  S'*)  +  |($^  ■S'*)rc 

=  {Q,S’')  +  K{HN,S'')rs  VS'*  GW'* 


(3.2) 


and 

Ka(fl\$'*)-hc(u,i?\$'*)-i(fl*,T'*-T„)r^  =  0  Vi^'*GW^  (3.3) 

Although  the  optimality  system  is  hnear,  the  coupling  of  u  and  $  in  the  two  equations  make  the 
derivation  of  error  estimates  nontrivial.  It  turns  out  to  be  convenient  to  apply  the  Brezzi-Rappaz-Raviart 
theory  (see  [2],  [5],  and  [7])  to  obtain  error  estimates. 

Theorem  3.1-  Let(T,^)  ond(r'*,$'*)  be  the  solutions  of  (2.9)-(2.10),  and  (3.2)-(3.3),  respectively.  Assume 
that  T,  $  G  n  for  some  I  <m<  k;  also  assume  that  (3.1)  holds.  Then, 

<  Cmax  1^,  l|  /l”*(||Q|lm-l  +  l|^f7v||r;v,m-l/2  +  l|Td||rc.m-l/2)  , 
where  C  is  independent  ofh,  6,  7,  T,  and  $.  ® 


A  simple  iterative  algorithm  for  solving  (2.10)-(2.11)  can  be  defined  as  follows: 
choose  ; 


for  n  =:  1, 2, . . . ,  solve  for  from 

Ka  (T<">,S)  +  c(u,T<"),S) 


=  +(<5.5)  +  «(/fAr,S)r^  VS  €  HhiQ) ; 

d 


(3.4) 


then  solve  for  from 


K  aiR,  *("))  +  c(u,  R,  $(“))  =  i(i2,  T<")  -  r<,)rc  V  iZ  €  ff],  (£i) . 
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(3.5) 


Of  course,  ultimately,  this  algorithm  has  to  be  carried  out  in  a  discretized  version,  such  as  one  using  a  finite 
element  method. 

The  convergence  of  this  algorithm  can  be  proved  as  a  result  of  the  observation  that  it  is  effectively  a 
gradient  method  for  the  following  minimization  problem:  find  g  €  L^iTc)  such  that  K(g)  :=  J(T{g),g)  is 
minimized  where  T{g)  6  ^h(^)  “  defined  as  the  solution  of  (2.1). 

Theorem  5.2-  lei  be  He  solution  <>/ (3.4)-(3.5)  and  (T,$)  the  solution  o/(2.9)-(2.10).  Then, 

Ji’*) -*T  in  and  $  in  .ff],(fl)  as  n-*oo.  I 

4.  Computational  Examples 

Let  Q  C  be  the  unit  square  (0, 1)  x  (0, 1).  Let  F  =  Fc  U  Fp  U  U  To  be  shown  as  in  Figure  1. 


Figure  1.  Computational  domain 


The  finite  element  spaces  are  chosen  to  be  piecewise  linear  elements  on  a  triangular  mesh.  All  the 
numerical  results  make  use  of  the  following  parameters  and  data: 


parameters: 


/c  =  l/0.73;  A5  =  1;  Wc^l;  w  =  l/2; 


boundary  data: 


r=l  onFp,  -3— =  0  onF^vUFo; 
on 


heat  source: 


Q  =  —  ^  ^9;r-  cos(3tx)  cos^(Ty)  —  sin*(^Tz)  cos(2Ty)^ 


+  (1  -  y*)-^  sin(3a-z)  cos=*(n-y) ; 
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velocity  profile: 


u  =  (1-J^,0), 


elastic  extra  stress  profile: 


pressure  profile: 

In  functional  (1.10),  we  choose 


Td  =  3.5 . 


For  the  data  given  above,  the  exact  solution  of  the  uncontroled  problem,  i.e.,  for 

^  =  0  onFc, 

an 

is  given  by  T  =  5  sin^(|T;r)  cos^(xy)  +  1. 

We  compare  the  temperature  distribution  in  the  uncontroled  case  with  the  optimal  temperature  distri¬ 
bution  in  the  controled  case  for  which 

dT  - 

^  =  ,  o.rc, 

where  g  is  the  control  such  that  (1.10)  is  minimized.  Approximations  to  the  optimal  state  and  co-state 
are  computed  from  (3.2)-(3.3);  the  approximate  optimal  control  is  then  obtained  from  (2.8),  t.c.,  g^  = 
—(1/5)$^ Ire-  All  of  the  computational  results  shown  below  were  obtained  with  the  use  of  a  mesh  size 
A  3=  Of  course,  calculations  with  varying  mesh  sizes  were  performed.  Since  these  merely  verified  the 
error  estimates,  we  do  not  report  on  them  here. 

Specifically,  Figures  2-4  deal  with  the  following  cases: 

1.  Exact,  uncontroled  temperature  and;  Optimal  temperature; 

2.  Exact,  uncontroled  temperature  and;  Optimal  temperature; 

2.  Optimal  boundary  control. 


Figure  2.  Temperature  contours:  uncontroled  and  controled. 
Optimal  boundary  control  on  Fc- 
(Fc  is  the  top  boundary  segment.) 

{To  is  the  right  boundary  segment.) 
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Figure  3.  Temperature  surfaces;  uncontroled  2uid  controled. 
(Fc  is  the  top-right  boundary  segment.) 

(Fo  is  the  lower-right  boundary  segment.) 
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Figure  4.  Optimal  boundary  control  on  Fc- 
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A  New  Condition  of  the  Exact  Boundary 
Controllability  for  Parabolic  Equations 


Yanzhao  Cao,  Max  Gunzburger  and  James  Turner 
April  26,  1997 


1  Introduction 

The  purpose  of  this  paper  is  to  investigate  the  exact  controllability  problems  for  systems  modeled 
by  linear  parabolic  differential  equations.  This  problem  has  been  studied  extensively  in  the  last 
30  years.  A  very  significant  early  contribution  is  due  to  Egorov  ([2],  [3]).  Another  early  paper 
is  due  to  Gal’chuk  ([9]).  Many  developments  in  the  controllability  theory  of  the  linear  parabolic 
equations  and  hyperbolic  equations  are  due  to  Fattorini  and  Russell  ([4],  [5],  [19],  [20]).  Using 
a  harmonic  analysis  method,  they  obtained  results  such  as  approximate  controllability  and 
conditions  on  y  for  the  exact  controllability  of  systems  governed  by  hyperbolic  and  parabolic 
equations.  We  refer  to  [7],  [8]  for  more  recent  development  by  Purskov  and  Imanuvilov. 

In  this  paper  we  attempt  to  study  the  exact  controllability  problem  by  examining  the  limit 
behavior  of  the  corresponding  optimal  control  problem.  Using  this  method  we  obtained  a  new 
condition  on  the  exact  controllability  of  the  system.  We  believe  that  this  condition  is  at  least 
comparible  condition  obtained  in  [18].  The  paper  is  organized  as  follows.  In  §2,  we  state  the 
problem  and  some  results  concerning  linear  parabolic  differential  equations.  In  §3  we  give  a 
constructive  proof  for  the  approximate  controllability  of  the  system  using  our  method.  Finally 
in  §4  we  give  a  condition  on  the  exact  controllability  of  the  system. 


2  Statement  of  the  Problem,  Notation  and  Preliminaries 

Let  fi  be  a  bounded  domain  in  For  a  fixed  T,  let  Q  =  [0,  T]  x  and  S  =  (0,T)  x  dCl.  On 
the  domain  let  A  be  the  second  order  elliptic  differential  operator 

^y  =  -Y. 

ij=l  *  ^ 

We  will  assume  that  c  >  0  on  fi  and  that  the  matrix  (aij(x))  is  symmetric  and  positive  definite. 
Consider  now  the  parabolic  initial  boundary  value  problem 

§+Ay  =  0,  inQ, 

g  =  tx,  onS,  (2.1) 

y(0)  =  yo,  in  fl  for  t  =  0. 

Hnre  g  =  j=iO,ijni{-^)  where  n*  is  the  i  —  th  component  of  the  outward  unit  normal  on 
5Q.  We  say  that  the  system  is  exactly  controllable  for  a  given  state  y  if  there  exists  u  €  L^(S) 
such  that  the  solution  y  of  problem  (2.1)  also  satisfies 

yiT,-)  =  y- 

Define  a  cost  functional 

Je{u)  =  (2/(^>  “)  “  yi^if  dx  +  I  ^  “^(i.  x)  dxdt , 

where  y  =  y{t,x;u)  satisfies  (2.1).  Uc{t,x)  is  said  to  be  an  optimal  control  if  it  satisfies 

Mue)  =  inf  ■  (2.2) 

u£L^{Q) 

For  a  Banach  space  X,  define 

L2(0,r;X)  =  {/  :  (O.T)  X;  T  ||/(t)l|^dt  <  oo}. 

Jo 

It  is  a  Hilbert  space  with  the  norm 

II/11l2(o,T;X)  =  (^  \\fi't)\\xdt)‘^  ■ 

Let  L^{Q,)  and  L^(0f2)  be  the  usual  function  spaces.  We  shall  use  the  notation 


and 


<(f)^'ip>=  / 

JdCl 

for  the  inner  products  and  the  associated  norms  will  be  denoted  by  ||  •  ||  and  1*1,  respectively. 
We  also  define 

W{0,T)  =  {/;/  e  €  L2(0,r;if-i(n))}, 

where  ^  is  taken  in  the  sense  of  distributions  on  (0,T)  with  values  in  W(0,r)  is  a 

Hilbert  space  with  the  norm 

ll/llw(o,T)  (ll/lli2(o,r;Hi(n))  ■ 

For  p,  q  >  0,  let  also 

HP’^iQ)  =  l2(0,  T;  D  HP{0,T;  L2(Q)). 

These  spaces  are  described  in  [15] ,  and  their  norms  are  defined  by 

II/IIhP'9((3)  =  (ll/lli2(0,r;H9(n))  +  ll/ll/fP(0,T;L2(f2)))^- 

The  spaces  HP’^CE)  with  associated  norms  j]  •  ||hp.9{E)  are  defined  similarly  by  replacing  SI  by  E 
above. 

We  recall  that  if  5  >  5  andp  =  then  there  is  a  constant  c  such  that  for  any  /  €  HP'^{Q), 

(E)  —  (Q)’ 

and  if  q  >  g  ^  and  t  e  [0,T]  then 

ll/WllHp-i(n)  -  (Q) 

and 

l5l//p-i(sn)  ^  ^11^11^9.! (s)- 

Define  the  bilinear  form  a:  H^{Q)  x  H^{Q,)  — >  i?  by 

Note  that  from  the  properties  of  the  operator  A  we  have  that  there  exist  constants  ci,  C2  >  0 
such  that 

\a{4>,i^)\  <  Cil|(?!)||i||^||i 
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and 

a{(j),(l>)  >  C2||(^||? 

for  all  (f>,  Ip  e 

We  now  define  the  operators  for  weak  solutions  of  elliptic  boundary  value  problems.  For 
/  €  let  Tf  e  denote  the  unique  solution  of  the  problem 

a{Tf,cP)  =  {f,cP)  forcPeHH^)- 

r  is  a  linear  operator  on  Also  since 

{Tf,(P)  =  a{Tf,T<P)  =  {f,T^) 

for  any  f,(p€  A^(0),  T  is  self  adjoint  on  L^{Q.).  Next  we  introduce  the  eigenvalue  problem 

A(p  —  A(p,  in  Q, 

^  =  0,  on  an. 

ou 

It  is  well  known  that  this  problem  has  a  system  of  eigenfunctions  forming  a  complete 

orthonormal  set  in  with  eigenvalues  0  <  Ai  <  A2  <  •  •  •  <  Aj  — oo  as  j  ^  oo.  Note  that 

if  we  let  Uj  =  A“^,  then 

T(j)j  =  ^j4>j  j  —  1?2, .... 

W  now  consider  the  weak  formulation  of  problem  (2.1) 

{^,cP)  +  a{y,<P)  =  <u,<p>,  foT<peH^{n), 

y(o,-)  =  vo¬ 

lt  is  well  known  that  if  yo  6  L^(Q)  and  u  G  then  (2.3)  has  a  unique  solution  y  in  W(0,r), 

and  there  is  a  constant  C,  independent  of  yo  and  u,  such  that 

l|y||w(o,r)  <  C'(||u||  +  ||w||l2(£;)). 

If  u  =  0,  then  the  exact  solution  of  (2.3)  can  be  represented  by 

OO 

j=l 

The  following  lemma  can  be  found  in  [21]. 

Lemma  2.1  Let  (j)  G  for  — 1  <  p  <  1.  Then 

OO 

is  an  equivalent  norm  in  i?P(fi). 
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3  Approximate  Controllability 


3.1  An  operator  representation  of  the  terminal  state.  In  the  following  lemma  we  state 
the  optimality  system  for  the  problem  (2.2).  It  will  play  a  central  role  in  obtaining  the  main 
results  of  this  section  and  the  next  section.  From  now  on,  we  suppress  the  explicit  dependence 
on  X,  e.g.,  we  will  use  u{t)  to  mean  u{t,x). 


Lemma  3.1  ([12])  The  optimality  system  of  problem  (2.1)  is  given  by 


° 

^  =  0 
Pe  =  Ve-y 
Pe  =  -eUe 


in  Q, 
on  S , 

in  n  for  t  =  r , 
on  E . 


(3.1) 


The  weak  formulation  of  the  optimality  system  can  be  written  as  follows. 

{^1^) — 

+  0,(Pei4>)  — 
y.(0)  = 

Pe{T)  = 

To  represent  the  terminal  state  ye{T)  by  an  operator,  we  consider,  for  a  given  z  e  L^(f^), 
the  system 

for0  6if^(O), 

+  o,(Pl->4>)  =  0)  for  €  i^^(r2), 

yi(0)  =  0, 

Pl(T)  =  2:. 

This  system  has  a  unique  solution  (yi(t),pi(t))  €  1^(0, T). 

Let  E{t)  denote  the  solution  operator  for  the  problem 

=  0  on  Q , 

=  0  on  S , 

=  2  in  f2  for  f  =  0 , 

w{t)  =  E{t)z . 


dw 

dw 

du 


W 


i.e., 


\<Pe,4>>,  ior  4)  e 
0,  for  4)  G  H\Q.), 
yo, 

y  -  ye{T). 
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We  now  define  an  operator  R:  — >  L^{Q.)  by 

Rz  =  2/1  (T). 


(3.4) 


Comparing  (3.2)  with  (3.3),  we  have  that 

^R{y{T)  -y)  =  -E{T)v  +  y{T)  (3.5) 

or 

(e  +  R){y{T)  -  y)  =  e{y  -  E{T)v)  (3.6) 

3.2.  Properties  of  the  operator  R. 

Lemma  3.2  R  €  L(i/^(f2),H^(f2)).  Therefore  R  is  a  compact  operator  in  L^(Q). 

Proof:  Consider  the  mapping 

z-*Pi{t,-)  =  E(T-t)z. 

We  prove  that  this  mapping  is  continuous  from  L^(f2)  into  By  standard  theory  for 

parabolic  differential  equations  ([15])  this  mapping  is  continuous  from  ff^(fi!)  into  If 

z  €  then 

Jo  Jo  ^ 

.  oo 

i=i 

<  C'lkll-i- 

Hence  the  mapping  is  continuous  from  H~^{Q)  into  L^(Q).  By  interpolation  ([1])  we  have  that 
Ris  &  continuous  mapping  from  L^(f2)  to  By  (3.3),  we  have  that  ^|e  =  Pi|e.  Thus 

llilzlli  =  ||yi(T)||i  <  C||yi||Hi,2(Q)  <  C||z||. 

The  proof  is  complete. 

Lemma  3.3  R  is  symmetric  and  semi-definite. 

Proof:  For  given  z,  zq  €  ^^(Q),  let  (y,p),  (2/0, Po)  be  the  corresponding  solutions  of  (3.3), 
respectively.  Then 

(^,Po)  +  a{y,Po)  =<  P,Po  >, 
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(-^,p)  +  o.{yo,p)  =<  Po,p  >, 


a(y,Po)  =  (2/’^Po), 


a{yo,p)  =  (yO’^P)- 


We  therefore  have  that 


{■^y{s),Po{s))  +  iy{s),  ^Po(s))  =  {■^yo{s),p{s))  +  {yo{s),  ^P(s)) 


^(yo(s),p(s))  =  ^{y{s),Po{s))- 

Integrating  the  above  equality  from  0  to  T,  noting  that  y(0)  =  yo(0)  =0,  we  have  that 


(y(r),po(T))  =  £  ^^{y{s),po{s))ds 

=  ^  £yo{s)p{s))ds 


=  (j/o(r),p(T)) 


(ilz,2o)  =  (^20,  2). 


Thus  R  is  self-adjoint.  We  also  have  that 


bP  =  i-^y,p)  +  aiy,p) 

=  i-Ty^p)  +  (.y^-j:P) 


■iy,p)- 


Integrating  the  above  equality,  we  have  that 


{Rz,z)^  /  |ppds>0. 

Jo 

So  R  is  semi-definite.  The  proof  is  complete. 

Lemma  3.4  Assume  that  S  is  analytic.  Then  Keri?  =  0. 

Proof:  Under  the  given  assumption  we  know  that  p  is  analytic  in  Q  and  on  2.  Assume  that 
Rz  =  0.  Then  from  the  proof  of  Lemma  3.3  we  have  that 


=  (i?2,2)=  [  \pfds. 
Jo 
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Thus  p  =  0  on  E.  Hence  the  data  corresponding  to  the  Cauchy  problem  is  zero  on  S  which 
implies  according  to  the  Cauchy-Kowaleska  theorem  (since  p  is  analytic)  that  p  =  0  on  Q.  In 
particular  z  =  p{T)  =  0.  The  proof  is  complete. 


3.3  Approximate  Controllability 

Theorem  3.5  Let  y  G  and  {ye,Ue)  =  (ye(t),'Ue(t))  be  the  solution  of  the  optimality 

stystem  .  Then 

Ill/e  (r)  -y  II L2(n)  =  0. 

Proof:  By  Lemma  3.2,  Lemma  3.3  and  Lemma  3.4,  and  the  Gilbert-Schmidt  theorem,  the 
operator  R  has  a  system  of  eigenfunctions  with  eigenvalues  Ai  >  A2  >  •  •  •  >  Aj  — »  0  as 

j  — ♦  00.  Moreover,  forms  an  orthonormal  basis  in  For  y  G  L^{Q)  let 

00 

y  -  E{T)v  =  '^yjej. 


Then 


Using  (3.6)  we  have  that 


Thus 


R{y  -  E{T)v)  =  Xjyjej. 

j=l 


y,{T)-y  =  Y 


^  €  +  A/ 

J  =  l 


00  ^2y2 


\\ye{T)-mhm  = 

i=i  '  ^ 


^  e2y2  00 

s  E  vl 

j=N+i 


For  any  given  6  >  0,  there  exists  N  such  that 


00  c 

E  vl  <  i 

j=N+l 


For  this  fixed  N  there  exists  0  >  0  such  that 


,^(e  +  A,)2  2 
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for  0  <  e  <  /?.  Thus 

\\ye{T)-y\\Ua)<S 

for  0  <  €  <  /?.  The  proof  is  complete. 

Definition  System  (2.1)  is  said  to  be  approximately  controllable  if,  for  any  given  y  and 
5  >  0,  there  exists  a  control  u{t)  and  a  function  y{t)  such  that  {y{t),u{t))  is  a  solution  of  (2.1) 
and 

The  following  theorem  is  an  immediate  consequence  of  Theorem  3.5. 

Theorem  3.6  System  (2.1)  is  approximately  controllable. 

Remark  Approximate  controllability  is  well  known  for  linear  parabolic  differential  equations. 
Our  proof  is  different  in  that  it  is  a  constructive  proof. 


4  A  Condition  for  Exact  Controllability 

In  this  section  we  give  a  condition  on  y  for  the  exact  controllability  of  system  (2.1).  The  main 
tool  is  the  operator  R  introduced  in  §2.2. 

Theorem  4.1  Let  R  be  the  operator  defined  in  ^2.2  and  ({ej}^i,  {Aj}^i)  be  the  eigensystem 
of  R.  Assume  that  R~^  is  the  inverse  of  operator  R.  If 

y  -  E{T)yo  ^  D{R-^) 

then  system  (2.1)  is  exactly  controllable. 


Hence 


°°  1 
oo  ifl  , 

i=l  ^ 

I 

Let  {ejk}^i  be  a  sequence  such  that  Cfc  ->  0  as  fc  oo.  Prom  the  optimality  system  (3.1)  and 
the  proof  of  Lemma  2.3  we  have  that 


l|y£fc(^)ll//l-2(Q) 


<  +  llyo||i/i(a)) 

<  <^(II^II^^.1(q)  +  II2/oIIhi(p)) 

<  “ — ^llL2(n)  +  l|yo||Hi(n)) 

<  C'illyollHi(n)  +  C2. 


Thus  there  exists  a  subsequence  of  {cfc},  still  denoted  as  {e^},  and  y  €  such  that 


lim  ye^  =  y  weakly  in 

fc-+CX> 

Using  the  trace  theorem  ([1])  and  the  compact  embedding  theorem  ([1])  in  Sobolev  spaces  we 
also  have  that  there  exists  u  €  L^(E)  such  that 

lim  =  u  in  L^(E). 

k—KX> 

Also  by  the  trace  theorem  ([1])  we  have  that  y{T)  =  y  and  y(0)  =  yo-  We  now  prove  that  {y,u) 
is  a  solution  of  (2.1).  Since  (ytfc(i),  «£*.(<))  is  a  solution  of  (2.2) 

{^,4,)  +  a{y,,,<l>)  =<  for  <j>  €  H\Q). 

Passing  to  the  limit  in  the  above  equality  we  obtain 

(§,</>) +a(y,.^)=<u,0>, 
at 

Thus  (y,u)  €  x  L^(S)  is  a  solution  of  (2.1).  This  proves  that  u  is  an  exact  control.  The 

proof  is  complete. 

Remark  4.1  Prom  the  proof  of  Theorem  4.1  we  see  that  if  (4.1)  holds,  then 


WyeiT)  —  y||L2(fi)  <  Ce. 


It  is  easy  to  show  that  (4.1)  is  actually  a  necessary  condition  for  the  above  inequality  to  hold. 
In  fact,  assume  that  the  above  inequality  holds.  Then  by  the  proof  of  Theorem  2.8  we  have  that 


Me,T)-nhi^)  =  E 


e'^y] 


<Ce^. 


Thus 


°° 


<C. 


Letting  e  — >  0,  we  obtain  (4.1). 

Remark  4.2.  In  [18],  Russell  gave  a  condition  on  exact  controllability  of  system  (2.1)  using 
the  harmonic  analysis  method.  The  condition  is 


y  -  EiT)v  e 


where  C  >  0  is  a  constant  and  A  is  the  elliptic  operator  that  appears  in  (2.1).  For  d  =  1,  this 
condition  implies  that  y  —  E{T)v  is  at  least  a  C°°  function.  At  this  moment  it  is  not  clear  if  our 
condition  is  weaker.  However,  if  the  control  is  distribution  we  can  show  that 


D{R-^)  =  D(A)  =  H^{n) 


which  is  much  weaker  than 
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